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Abstract 



We show that the neutron and proton transition densities predicted by recent 
Quantum Monte Carlo calculations for A = 6, 7 nuclei are consistent with pion 
scattering from 6 Li and Li at energies near the A resonance. This has pro- 
vided a microscopic understanding of the enhancement factors for quadrupole 
excitations, which were needed to describe pion inelastic scattering within the 
nuclear shell model of Cohen and Kurath. 

PACS numbers: 21.60.Ka, 25.80.Ek 
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Quantum Monte Carlo (QMC) methods have been successfully developed to predict the 
properties of low-lying states of light nuclei starting with realistic two- and three-nucleon 
potentials [I]. While the reproduction of the energy spectra is the essential first step in 
such an effort, the dynamical content of the resulting nuclear wave functions must be tested 
against various reactions. This has been achieved so far mainly by considering electroweak 
processes [@||. In this paper we report on an additional test by using pion elastic and 
inelastic scattering. 

Let us first briefly review the status of our understanding of pion inelastic scattering at 
medium energies( 80 MeV < < 300 MeV). Because of the excitation of the A resonance, 
the pion-nucleus interactions in this energy region are dominated by the strong absorption 
mechanism. Consequently, the pion-nucleus inelastic scattering leading to discrete final 
nuclear states can be described by the Distorted Wave Impulse Approximation (DWIA). 
This has been well established [{4)^1311 in very extensive investigations of the data from meson 
factories. Following the momentum-space approach [0], the inelastic scattering amplitude 
can be written as 

Tfifa, k ) =Jdk'J dk X ^ f (k')U fl {k\ k) X ^\(k) , (1) 

where the fc's are pion-nucleus relative momenta in the pion-nucleus center of mass frame, 
and the distorted waves x are generated from an optical potential which is adjusted to 
fit the pion-nucleus elastic scattering. 

The nuclear excitations are contained in the transition potential Up. It can be calculated 
from the ttN scattering t-matrix and nuclear transition form factors. The details are given 
in Ref. [0. To simplify the presentation, the spin-isospin variables will be suppressed here. 
Then the transition potential can be written as 

U fi (k', k) = t nn {k', k, u )fj^ ,2 {q) + t np (k', k, u )f){ 2 (q) , (2) 

where q = k' — k, t nn (t np ) is an appropriately parameterized pion-neutron (pion-proton) 
scattering amplitude, and lvq is the collision energy calculated from using the fixed-scatterer 
approximation. With = —1/2, 1/2 denoting neutron and proton respectively, the nuclear 
transition form factors are defined by 

f%{q)= Jdre-^p%(f), (3) 

with the transition densities defined by 

=< */ I \ E 6(r~ rl) ^ 2 ^^ I *< > , (4) 

i=l,A " 

where r 2 is the z-component of the nucleon isospin operator, and and ^ / are the initial 
and final nuclear states respectively. The transition dynamics can be better understood 
from the multipole expansion of transition densities 

P%(fi = E i W) < J f M f I JiKMiM > Fg %(r) , (5) 

KM ' 
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Here we recall the notation of Ref. J7j to include possible spin transitions S = 0, 1 in 
F^g t (r), where K denotes the orbital angular momentum transition and J = K + S is the 
total angular momentum transfer. 

In the DWIA study for lp-shell nuclei, the transition densities were calculated from 



the shell model of Cohen and Kurath [TJJ 



F KS,t 3 (r) =J2 A J(KS)ts(a,& f^ Ja ) 1/2 < (U/2)j Q || [Y K {r) x a s ]j \\ (l p l/2)j p > 

a/3 

x Rn a i a (r)Rn ,i (r) (6) 

with cr = 1, (J\ = a, and 

A JiKS )t s (a, P] fi) =< */ || [bit, x hl t3 \ J[KS] \\ % > , (7) 

where a(n a ,l a , j a ) denotes the single particle orbitals, R a {r) is the radial wave function, 
and and hg are the creation operators for the particle and hole states respectively. 

It was found that the pion inelastic scattering from lp-shell nuclei can be described 
by using the above shell-model input only when the quadrupole excitation component 
J(KS) = 2(20) is enhanced by a factor E N ~ 2. In Ref. 0, these enhancement factors 
were estimated from a systematic analysis of B(E2) transitions from lp-shell nuclei. For 
proton excitation, this enhancement factor is consistent with what is needed for explaining 
B(E2) values. For Z = N nuclei, one can assume that the neutron excitation also has the 
same enhancement because of isospin invariance. However, for the N ^ Z nuclei, such as 
7 Li, the enhancement factors for neutron excitations can not be obtained without making 
some additional assumptions. The predicted pion inelastic cross sections thus are not well 
justified theoretically. Furthermore, it would be desirable if the calculations do not include 
any enhancement factors. This however is very difficult, if not impossible, in practice within 
the shell model since the collective quadrupole excitations can only be described by a very 
large model space. 

In this work, we calculate the transition densities for a given multipolarity by using wave 
functions from the recent QMC calculations for light nuclei. The input, Eq.(6), to the DWIA 
calculations for pion inelastic scattering is then defined by the following matrix element 









<J f M f 


./,;./ .\/, .U > 



The QMC calculations use a realistic Hamiltonian containing the Argonne Vx% two- 
nucleon [ITSI and Urbana IX [HJ three-nucleon potentials, which we refer to as the AV18/UIX 
model. Both variational (VMC) and Green's function (GFMC) Monte Carlo calculations 
have been made for light nuclei Q. The AV18/UIX model reproduces the experimental 
binding energies and charge radii of 3 H, 3 He, and 4 He, in the numerically exact GFMC 
calculations, but underbinds 6 Li and 7 Li by 2-5%. The variational Monte Carlo (VMC) 
energies are 2% above the GFMC results for A = 3, 4 nuclei and 10% above for A = 6, 7. 
However, the known excitation spectra are well reproduced by both the VMC and GFMC 
calculations, as are the charge radii. The VMC and GFMC calculations also produce very 
similar one-body densities, while two-nucleon density distributions differ by less than 10%. 
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The VMC wave functions have been used successfully to describe the elastic and tran- 
sition electromagnetic form factors for 6 Li Q] without introducing effective charges. They 
have also given an excellent absolute prediction for the spectroscopic factors in 7 Li(e, e'p) 
reaction ||. Consequently we expect the VMC wave functions to give a good estimate for 
both the elastic and transition densities required in pion scattering calculations. 

The variational wave function for A = 6, 7 nuclei used here is the trial wave function, 
that serves as the starting point for the GFMC calculations. It has the general form 



|*T> 



E uW 

i<j<k 



SU(l + U i: 



i<j 



(9) 



where £/y and are two- and three-body correlation operators and the Jastrow wave 

function is given by 



\*j)=m n fk n /. 



' ijk 

i<j<k<4 «<i<4 



n fs P ( r ki) 

k<A<KA 



E 

LS 



LS[n] 1234:56. 
4<l<m<A 




(10) 



The S and A are symmetrization and antisymmetrization operators, respectively. The cen- 
tral pair and triplet correlations fxyfrij) and f^ k are functions of relative positions only; the 
subscripts xy denote whether the particles are in the s- or p-shell. The \^A(LS[n]JMTT 3 )) 
is a single-particle wave function with orbital angular momentum L, spin S, and spatial 
symmetry [n] coupled to total angular momentum J, projection M, isospin T, and charge 
state T 3 : 



\$ A (LS[n]JMTTi 



3jl234:56...A/ = i$ Q (oooo) 1234 n 

4</<A 



ad i 



IJ Y lm (n al )} LML[n] x [ JJ Xi{\m s )}sM s 



n m¥s)]tt 3 ) ■ 



(ii) 



4<1<A 



4<1<A 



JM 



4</<A 



Particles 1-4 are placed in an a core with only spin-isospin degrees of freedom, denoted by 
$ a (0000), while particles 5-A are placed in p-wave orbitals 4>p S {Rai) that are functions of the 
distance between the center of mass of the a core and particle I. Different amplitudes @LS[ri\ 
are mixed to obtain an optimal wave function by means of a small-basis diagonalization. 
For 6 Li, the (J n ;T) = (l + ;0) ground state is predominantly a 3 S[2] amplitude, with small 
admixtures of 3 D[3] and 1 P[11] components, while the (3 + ; 0) first excited state is pure 3 D[3]. 
For 7 Li, the (J 7 "; T) — (| ; |) ground and (\ ; \) first excited states are predominantly 2 P[3], 
with small admixtures of 2 ' 4 P[21], 2 - 4 D[21], and 2 S[111] components. The (J^;T) = (| _ ; |) 
and (| ; |) excited states are predominantly 2 F[3], again with small admixtures of 2,4 P[21] 
and 2 ' 4 D[21] components. Mixing parameter values are given in Ref. 
The two-body correlation operator is defined as: 



U, 



E 

p=2,6 



II fijk( r ik,r jk ) 



U, 



ij ) 



(12) 



4 



where the Of~ 2,6 = Tj • Tj, cr^ ■ crj, cr^ • (TjT{ ■ Tj, SV,-, and SijTi ■ Tj. The six radial functions 
f ss (r) and u p=2 fi(r) are obtained from two-body Euler-Lagrange equations with variational 
parameters [|T7j]. The f sp and fp p s ™ correlations are similar to / ss for small separations, but 
include parameterized long-range tails. The parameters used in constructing these two-body 
correlations, as well as the description of the three-body correlation operator U™ 1 and the 
operator- independent three-body correlations f^ k and ffj k are given in Ref. 

In Ref. [0 it was found that inelastic transitions induced by pion scattering are dominated 
by the quadrupole transition J(KS) = 2(20). In our QMC calculations, we evaluate the 
quadrupole transition density: 



l2J f + l(* JfMf | £ i=M J(r - rdrfYj^^f^ \ * JiMi ) 
pE ^ r) ~ (JfMf | Ji2MiM) • {U) 

These neutron and proton transition densities are shown in Fig. |l] for four transitions in 6 Li 
and 7 Li. The integrated B(E2 1,ts) values, 

B(E2 th) ^^*^ (14) 

are given in Table p], where they are compared to the experimental proton values obtained 
from (e, e') scattering and Coulomb excitation experiments |l8|-p0"||. Our evaluations are 
made with 160 000 Monte Carlo samples for transitions in 6 Li and 120 000 Monte Carlo 
samples for transitions in 7 Li. This number of samples is sufficient to give statistical un- 
certainties that are as small or smaller than the errors on the experimental B(E2) values. 
Other J(KS) transition amplitudes can contribute to pion inelastic scattering, particularly 
in the 6 Li case, and are evaluated in a similar manner. 

We see from Fig. [I] and Table | that the predicted differences between the neutron and 
proton excitations in 7 Li are very significant. Such differences can be most effectively verified 
by using an important characteristic of pion scattering at energies near the A excitation. At 
a typical energy E w = 164 MeV, one finds that the irN amplitude in Eq.(2) has an interesting 
ratio | t^+p/t^+n | = | t^-n/t^-p |~ 3. Consequently, the 7r+ scattering is dominated by the 
proton excitations while n~ scattering is dominated by neutron excitations. The agreement 
with both the 7r + and 7r _ data will be a nontrivial test of the QMC wave functions. Thus the 
present study is complementary to that of Ref. || using electron scattering which mainly 
probes the proton excitations. 

We first investigate pion scattering from 6 Li. Here data for E n = 100, 180 and 240 
MeV [[L3| are available for testing the energy-dependence of our predictions. The pion optical 
potential and the irN t-matrix are taken from Ref. 0. The nN amplitudes we employ are 
taken from the Karlsruhe-Helsinki analysis [EI], and differ only slightly with the more recent 



VPI analysis [^], mainly in the Sn partial wave (as discussed in Ref. @]). This partial 
wave and the other non-P 33 partial waves are much weaker than the P33 channel in the 
energy region of interest near the A excitation. For the present exploratory investigation, 
we therefore do not make any effort to improve the optical potential employed in Ref. [0. 
Such an improvement is probably needed in the future when the data at low energies and 
for the spin observales are investigated. 

With transition densities calculated from QMC, there are no adjustable parameters in 
our DWIA calculations. In Table |, we see that the calculated B(E2) for the transition to the 
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(3 + , T = 0) state is in excellent agreement with the data. This is a significant improvement 
over the shell-model prediction which required a large enhancement factor E p = E n = 2.5 to 
reproduce the B(E2) data, as discussed in Ref. f7|. We thus expect a similar improvement 
in pion scattering calculations. 

Our results for 6 Li are displayed in Fig. |2|. We see good agreement with the differential 
cross sections for elastic scattering. However, the discrepancies in reproducing the diffrac- 
tive minima at 180 MeV indicate some deficiencies of the simple optical potential we have 
employed. We also see general agreement with the inelastic scattering to the (3 + ;T = 0) 
excited state, although some noticeable discrepancies are seen, particularly at 100 MeV. 
Nevertheless, it is fair to say that our results agree with the differential cross sections to a 
very large extent in both absolute magnitude and energy-dependence. The agreement seen 
in Fig. Q is consistent with the B(E2) values listed in Table |. The overall agreement is 
not surprising in view of the ability of the QMC wave functions to reproduce the elastic 
and transition form factors in electron scattering experiments 0. To further improve the 
agreement with the data and to account for the spin observables, it would be necessary to 
improve the reaction model. For example, we may have to consider coupled- channel effects 
and refine the optical potential, but this is beyond the scope of the present investigation. 

We next investigate the very old data for 164 MeV 7r + and tt~ inelastic scattering from 
a 7 Li target |24|]. The final states we consider are (J n ;T) = (| ; |) at 0.478 MeV, (| ; |) 
at 4.63 MeV, and (§"; ~) at 6.68 MeV. Our results are shown in Fig. |3|. We see that the 
predicted cross sections (solid curves) are in excellent agreement with the data. In the 
same figure, we also show the contributions obtained using just the proton excitations. 
The agreement with both the n + and 7r~ data is evidently due to the delicate interplay 
between the neutron and proton excitations. If the shell model input given in Ref. |7j] is 
used without the enhancement factors E n = 1.75 and E p = 2.5 for the quadrapole transition 
J(KS) = 2(20), the predicted cross sections will be a factor of about 5 lower than the data 
for all of the cases considered here. 

In conclusion, we have performed calculations of pion scattering from 6 Li and 7 Li using 
the nuclear transition densities predicted by the recent QMC calculations for light nuclei 
starting with realistic two-nucleon and three-nucleon potentials. The predicted cross sections 
are in very good agreement with the data. In contrast with the previous calculations using 
densities from nuclear shell model, the calculation does not include any enhancement factors. 
Because of the strong isospin dependence of the ttN scattering t-matrix, the present inves- 
tigation has probed critically the predicted neutron transition densities which are not well 
tested in electron scattering studies. Our results suggest that the wave functions predicted 
by the QMC calculations are accurate for investigating various nuclear reactions. 

It is highly desirable to extend the present work to re-investigate the very extensive data 
of pion-nucleus scattering on larger lp-shell targets and other more complex processes such as 
pion absorption and double-charge-exchange reactions. We expect that QMC wave functions 
for A = 9, 10 nuclei will become available in the next year. With the nuclear correlations 
correctly accounted for by using the wave functions predicted by QMC calculations, one 
now can hope to resolve many long-standing problems in intermediate-energy pion-nucleus 
reactions. 
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TABLES 



TABLE I. B(E2 ^,t$) values in e 2 -fm for different transitions in 6 Li and 7 Li. Experimental 
values are from Refs. |l8| , p!s[| . Experimental uncertainties and Monte Carlo sampling errors are 
given. 



A Z(J, - 




Experiment 
P 


P 


QMC 


n 


6 Li(0+ - 


-3+) 


21.8±4.8 


21.1±0.4 




21.1±0.4 


T Li(§~ - 


-n 


7.59±0.10 


5.7±0.1 




16.5±0.3 


7 Li(f- 


-n 


15.5±0.8 


13.2±0.2 




34.6±0.5 


7 Li(f- - 


-n 


4.1±2.0 


2.3±0.1 




5.4±0.2 
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FIGURES 

FIG. 1. The transition densities, p^'f, for 6 Li and Li. 

FIG. 2. Differential cross sections for 6 Li(7r,7r) and 6 Li(-7r, tt') scattering at multiple pion ener- 
gies. The data are taken from Ref. fi~3|| . 

FIG. 3. Differential cross sections for 7 Li(-7r, tt') scattering at E n = 164 MeV. The data are 



taken from Ref. 24] 
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Fig. 1 Lee & Wiringa 
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Fig. 2 Lee & Wiringa 
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Fig. 3 Lee & Wiringa 
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